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ABSTRACT 

Future galaxy redshift surveys aim to measure cosmological quantities from the galaxy power 
spectrum. A prime example is the detection of baryonic acoustic oscillations (BAOs), provid- 
ing a standard ruler to measure the dark energy equation of state, w(z), to high precision. The 
strongest practical limitation for these experiments is how quickly accurate redshifts can be 
measured for sufficient galaxies to map the large-scale structure. A promising strategy is to 
target emission-line (i.e. star-forming) galaxies at high-redshift (z ~ 0.5-2); not only is the 
space density of this population increasing out to z ~ 2, but also emission-lines provide an 
efficient method of redshift determination. Motivated by the prospect of future dark energy 
surveys targeting Ha emitters at near-infrared wavelengths (i.e. z > 0.5), we use the latest 
empirical data to model the evolution of the Ha luminosity function out to z ~ 2, and thus 
provide predictions for the abundance of Ha emitters for practical limiting fluxes. We caution 
that the estimates presented in this work must be tempered by an efficiency factor, e, giving 
the redshift success rate from these potential targets. For a range of practical efficiencies and 
limiting fluxes, we provide an estimate of uPq,2, where n is the 3D galaxy number density 
and Prj.2 is the galaxy power spectrum evaluated at k — 0.2/iMpc~ . Ideal surveys must 
provide nPo.2 > 1 in order to balance shot-noise and cosmic variance errors. We show that a 
realistic emission-line survey (e = 0.5) could achieve nPo.2 = 1 out to z ~ 1.5 with a lim- 
iting flux of 10~ 16 erg s _1 cm~ 2 . If the limiting flux is a factor 5 brighter, then this goal can 
only be achieved out to z ~ 0.5, highlighting the importance of survey depth and efficiency 
in cosmological redshift surveys. 
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1 INTRODUCTION 

One of the greatest challenges the current generation of cosmolo- 
gists faces is to understand the physics underlying the apparent ac- 
celeration of the expansion of the Universe (e.g. Riess et al. 1998; 
Peiimutter et al. 1999). Contemporary models favour the influence 
of a dark energy that has come to dominate the energy density of 
the universe during the last 8 billion years. Unfortunately dark en- 
ergy is outside the realm of the standard model, and requires new 
physics to explain. Nevertheless, many mechanisms have been pro- 
posed, and the potential for establishing which (if any) is correct 
experimentally, has caused great fervour amongst the astronomical 
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community over the past decade. The reward for investing a large 
amount of effort into determining the physics of dark energy is of 
course a profound advancement of our understanding of the funda- 
mental nature of the universe. 

A range of dark energy models exist (see Peebles & Ratra 
2003 and Copeland et al. 2009 for reviews), however the two most 
prominent scenarios attribute the accelerating expansion to (a) a 
'cosmological constant' (A) analogous to a non-zero quantum me- 
chanical vacuum energy that has now come to dominate the over- 
all energy density of the universe (but 120 orders of magnitude 
smaller than the value predicted by quantum physics); or (b) a dy- 
namic scalar field ('quintessence') which varies with both time and 
space. Both models require general relativity to hold on cosmo- 
logical scales. A third alternative to explain the acceleration is the 
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Table 1. Parameters of the luminosity functions used to derive the empirical model of Ha counts. All Schechter function parameters have been corrected to a 
common fiducial cosmology (Ho = 70kms~ 1 Mpc , n m = 0.3, Ha = 0.7). 



Reference 


2 
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logL*(erg s _1 ) log</)*(Mpc~ 3 ) 


ct 


EW (A) 
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Gallego et al. (1995) 


<0.045 


41.87 


-2.78 


— 1.3 


>10 


UCM survey 


Shioya et al. (2008) 


0.24 


41.94 


-2.65 


— 1.35 
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Narrowband 0.815/^m 


Yan et al. (1999) 


1.3 ± 0.5 
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-1.35 
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Geach et al. (2008) 


2.23 ± 0.03 
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-2.84 


-1.35 


>12 


Narrowband 2.121/xm 



failure of general relativity on large scales, such that the gravity 
theory itself needs to be modified (e.g. Dvali et al. 2000). 

One way of distinguishing A from quintessence is to measure 
the evolution of the expansion of the universe, which is controlled 
by the dark energy equation of state, w(z); the ratio between dark 
energy pressure P and density, p. For A models, p = — P/c 2 for 
all time, such that w(z) — — 1. Detecting a varying w(z) would be 
a possible indication for quintessence. However, if gravity is modi- 
fied, such behaviour could be just an indirect effect of the failure of 
general relativity. Such degeneracy between dark energy and mod- 
ified gravity can be lifted only by measuring the growth rate of 
cosmic structure f(z) (or its integral G(z)), which is governed by 
the interplay between the strength of gravity and the expansion rate 
of the Universe. Thus, measuring w(z) and f(z) as a function of 
redshift represents the most powerful combination of observational 
probes to distinguish among competing models (Guzzo et al. 2008, 
Wang et al. 2008). This is the primary goal of current and future 
dark energy surveys (Albrecht et al. 2009). 

Although the accelerated expansion was discovered using ob- 
servations of Type la supernovae (Riess et al. 1998; Perlmutter et 
al. 1999), results from these observations are now dominated by 
systematic errors (e.g. Hicken et al. 2009). Future studies of dark 
energy therefore aim to exploit different observations and of those 
proposed, the use of Baryon Acoustic Oscillations (BAO) as stan- 
dard rulers appears to have the lowest level of systematic uncer- 
tainty (Albrect et al. 2006). BAO are a series of peaks and troughs 
in the power spectrum, which quantifies the clustering strength 
of matter as a function of scale. They occur because primordial 
cosmological perturbations excite sound waves in the relativistic 
plasma of the early universe: when the plasma breaks down at 
recombination, the radiation can be observed as the Cosmic Mi- 
crowave Background (CMB), while the fluctuations in the baryonic 
material give rise to BAO (Silk 1968, Peebles & Yu 1970, Sun- 
yaev & Zel'dovich 1970, Bond & Efstathiou 1984, 1987, Holtzman 
1989). The BAO signal is on large-scales, which are predominantly 
in the linear regime today. It is therefore expected that BAO should 
also be seen in the galaxy distribution (Goldberg & Strauss 1998, 
Meiksin, White & Peacock 1999, Springel et al. 2005, Seo & Eisen- 
stein 2005, White 2005, Eisenstein, Seo & White 2007, Kazin et al. 
2009), and can be used as a standard ruler, leading to measurements 
of the angular diameter distance Da(z) and the Hubble expansion 
rate H(z), and therefore w(z) (Seo & Eisenstein 2003, Blake & 
Glazebrook 2003, Hu & Haiman 2003, Wang 2006). 

The acoustic signature has now been convincingly detected at 
low redshift (Percival et al. 2001, Cole et al. 2005, Eisenstein et 
al. 2005, Huetsi 2006) using the 2dF Galaxy Redshift Survey (2dF- 
GRS; Colless et al. 2003) and the Sloan Digital Sky Survey (SDSS; 
York et al. 2000). Further analyses of the SDSS have led to com- 
petitive constraints on cosmological models (Percival et al. 2007, 
Gaztanaga et al. 2009, Percival et al. 2009). Ongoing spectroscopic 



surveys aiming to use BAO to analyse dark energy include the 
Baryon Oscillation Spectroscopic Survey (BOSS; Schlegel, White 
& Eisenstein 2009), the Hobby-Eberly Dark Energy Experiment 
(HETDEX; Hill et al. 2008) and the WiggleZ survey (Glazebrook 
et al. 2007). Ongoing photometric surveys such as the Dark Energy 
Survey (DES: [http : / /www . darkenergysurvey ■ or g), the 
Panoramic Survey Telescope & Rapid Response System (Pan- 
STARRS: |http : / /pan-starrs ■ if a ■ hawaii . edu[ l aim to 
find BAO using photometric redshifts. 

The power spectrum (or correlation function) of the galaxy 
distribution also contains key information on the growth rate of 
structure f(z) (Kaiser 1987). This produces large-scale motions 
towards density maxima, that contribute a peculiar velocity com- 
ponent to the measured galaxy redshifts used to reconstruct cosmic 
structure in 3D. The net effect is to produce an anisotropy in the 
power spectrum that can be measured to extract an estimate of the 
growth rate f(z), modulo the bias factor of the galaxies being ob- 
served. The importance of this well-known effect in the context of 
dark energy has become evident only in recent times, when redshift 
surveys of sufficient size at z ~ 1 have started to become available 
(Guzzo et al. 2008). Thus, a redshift survey of galaxies provides us 
with the ability to obtain an estimate of both key probes of cosmic 
acceleration, the expansion rate and the growth rate. 

In order to reduce shot-noise and cosmic variance in 'pre- 
cision' measurements of BAO and redshift distortions, the ul- 
timate observational challenge is to accurately measure a large 
number (tens or hundreds of millions) of redshifts for galax- 
ies spread over a significant interval of cosmic time, spanning 
the transition from matter domination to dark energy domina- 
tion in the universe, and covering the majority of the extragalac- 
tic (|b| > 20°) sky, ~2 x 10 4 square degrees. Such a sur- 
vey can be conducted using a dedicated survey telescopes from 
a space platform, as proposed by the Joint Dark Energy Mis- 
sion (JDEM: |http : / / jdem . gsf c . nasa ■ gov| > and European 
Space Agency's Euclid and SPACE satellite mission concepts 
( |http ://sci.esa. int/euclid ; Cimatti et al. 2009). Some 
of the ongoing and planned future BAO surveys, such as Wig- 
gleZ, will target emission-line galaxies - i.e. generally star-forming 
galaxies with easily identifiable redshifts. The goal of this work is 
to make a prediction for the abundance of Ha emitting galaxies that 
these dark energy surveys can expect using the existing empirical 
evidence of past and recent Ha surveys out to z ~ 2. 

In this work we use empirical data to build a simple phe- 
nomenological model of the evolution of the Ha luminosity func- 
tion (LF) since z ~ 2, and therefore predict the number counts 
of Ha emitters in redshift ranges pertinent to future dark energy 
surveys (the empirical model can also be used as a fiducial point 
for semi-analytic predictions for the abundance of star forming 
galaxies, e.g. Baugh et al. 2005; Bower et al. 2006; Orsi et al. 
2009 in prep). In Section 2 we describe the model, list the prin- 
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Figure 1. Evolution of the Ha luminosity function, assuming our simple 
model of L* oc (1 + z)Q out to z = 1.3 and no evolution to z < 2.2. 
The panels show the LF at z ~ 0, 0.2, 0.9 & 2.2, with observational data 
overlaid (all data has been corrected to the same fiducial cosmology used 
throughout this work and not corrected for extinction). Note that not all 
of the observational data shown here was used to construct the model (see 
§2.1), however the model is a good representation of the observed LFs out 
to z ~ 2. The largest discrepancy occurs at z ~ 1, where there is some 
scatter between different surveys. However, in part, this is due to the mixture 
of survey strategies, and cosmic variance in the small fields observed. The 
model LFs have been truncated at the luminosity limit corresponding to a 
flux of 10 _1B erg s — 1 cm -2 at each epoch (vertical dotted lines). Although 
there are hints that the faint-end slope is steepening out to z ~ 1, in the flux 
regime of practical interest this does not have a significant impact on our 
counts (also see §2.2.1 & Figure 3). 

cipal predictions and draw the reader's attention to some important 
caveats. In Section 3 we discuss the implications of the number 
count predictions on planned dark energy surveys, and in Section 4 
we comment on the relevance of cosmological surveys in the near- 
IR from a terrestrial base. For luminosity estimates, throughout we 
assume a fiducial cosmological model of Ho = 70 km s _1 Mpc^ 1 , 
n m = 0.3 and ^ A = 0.7. 



Figure 2. A comparison of the predicted number counts of Ho emitters 
from the simple model to observed counts integrated over the redshift range 
0.75 < z < 1.90. The shaded region indicates the \a uncertainty on 
the model counts. We compare to the observational data of the (slitless 
spectroscopic) surveys of McCarthy et al. (1999), Hopkins et al. (2000) 
and Shim et al. (2009), where the integrated counts have been calculated 
from the respective luminosity functions, uncorrected for dust extinction 
(so the counts include incompleteness corrections specific to each survey). 
Note that all-sky redshift surveys are unlikely to probe below flux limits of 
~10 — 16 erg s _1 cm -2 , where uncertainties due to the poorly constrained 
faint-end slope become more important to the count predictions (see §2.2.1 
for more details). 

(e.g. Gallego et al. 1995; Ly et al. 2007; Shioya et al. 2008). This 
is a favourable line to target at high-redshift because it is the least 
affected by extinction (compared to, say, [Oil]). The shape of the 
Ha LF in the local Universe is well characterised, and over the past 
decade, near-infrared surveys have tracked the evolution of the LF 
out to z ~ 2 (McCarthy et al. 1999, Yan et al. 1999, Hopkins et al. 
2000; Moorwood et al. 2000). Furthermore, the increasing feasibil- 
ity of statistically significant wide-field Ha surveys at high-redshift 
have vastly improved our picture of how the Ha luminosity func- 
tion has evolved over the past 8 Gyr (e.g. Geach et al. 2008; Shim 
et al. 2009; Sobral et al. 2009). 



2 A SIMPLE MODEL OF THE EVOLUTION OF THE Ha 
LUMINOSITY DENSITY 

Fortuitously for dark energy surveys, the global volume averaged 
star formation rate increases steeply out to z ~ 2, and flattens 
(or perhaps gently declines) towards earlier epochs (e.g. Lilly et 
al. 1995, Hopkins 2004). This will work in favour of dark energy 
surveys, provided the shape of the LF is reasonably well under- 
stood. Locally, star forming galaxies can be easily selected using 
the well calibrated and 'robust' Ha emission line at A = 6563A 



2.1 Empirical fit 

Throughout this work, we assume the conventional form of the lu- 
minosity function holds at all epochs - the Schechter function: 

4>(L)dL = cj>*(L/L*) a exp(-L/L*)d(L/L*) (1) 

In Table 1 we list the Schechter function parameters derived from 
four Ha surveys spanning < z < 2, chosen for their similarity 
in fitting (all find or fix the faint-end slope a = —1.35) and equiv- 
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Figure 3. Predicted redshift distribution d N/dz of Ha emitters for limiting 
fluxes of 1-5 X 10~ 16 erg s~ 1 cm~ 2 (thick to thin lines). Note that the tran- 
sition between L* evolution and non-evolution at z = 1.3 introduces the 
sharp fall-off in counts towards high-z. For comparison, we also show the 
redshift distribution for the same L* evolution and fixed </>*, but allowing 
the faint end slope to steepen monotonically from —1.35 at z = to —1.6 
at z = 2. The impact this change has on the predicted counts in the flux 
limits of practical interest is negligible, and (as expected) more pronounced 
at fainter limits. 



alent width cuts, generally EWo > 10A. Note that there is very 
little evolution in the LF between z ~ 2 and z = 1.3 (Yan et al. 
1999; Geach et al. 2008), although both of these surveys assume 
a fixed faint-end slope of —1.35 similar to that found in the lo- 
cal Universe (necessitated by the depths of these surveys). In com- 
parison, by z ~ 0, the characteristic luminosity L* has dropped 
by an order of magnitude (Gallego et al. 1995). The evolution of 
the space density normalisation <f>* is harder to model - the val- 
ues listed in Table 1 imply little evolution (compared to L*) with 
((f)*) — 1.7 x 10~ 3 Mpc~ 3 . However, other surveys have derived 
a larger range of <f>* (e.g. Sobral et al. 2009), probably in part due 
to cosmic variance effects, and the inherent degeneracy in LF pa- 
rameter fitting. The latter is the main reason we chose surveys with 
very similar fitting techniques; an attempt to mitigate the impact of 
different survey strategies on our model. 

With this in mind, the model presented here assumes evolution 
only in L* , and the faint end slope is held fixed at a = —1.35 (we 
assess the impact of this assumption in §2.2.1). Given the strong 
luminosity evolution out to at least z = 1.3, and weak evolution 
beyond to z ~ 2, we model the L* evolution as (1 + z) Q over 
< z < 1.3 (the median redshift of the HST/NICMOS grism 
survey of McCarthy et al. [1999]). At z > 1.3 we freeze evolution, 
and assume this is valid out to the limit of current Ha observations 
(z = 2.23). The best fit L* evolution is then derived as: 



We estimate the uncertainty in Q via a bootstrap-type simulation; 
re-evaluating the fit 10,000 times after re-sampling each L* in Ta- 
ble 1 from a Gaussian distribution of widths set by the L* la un- 
certainty. Note that L* has not been corrected for intrinsic dust 
extinction (a canonical Ana = 1 mag correction is generally ap- 
plied when deriving star formation rates, although this could in- 
crease at high luminosity). The luminosities have been corrected 
for [N II] contribution, typically of order ~30% (e.g. Kennicutt 
& Kent 1983). Note that this could be a conservative correction 
if there is a significant contamination from active galactic nuclei 
(AGN). With this in mind, Ha redshift surveys should aim for a 
spectral resolution that can resolve Ha/[N II]. Not only does this 
have a significant practical benefit, in that it aids redshift identi- 
fication, but also the secondary science impact of a large sample 
of Ha/[N II] ratios, and thus AGN selection would be extremely 
valuable. 

As described above, the choice of normalisation of the model 
is a source of uncertainty in the predicted counts. Since this pa- 
per is focused on predictions for dark energy surveys, which will 
target Ha emitters at z ~ 1, here we have taken the normalisa- 
tion of the model to be the average <j>* of the surveys of Yan et 
al. (1999), Hopkins et al. (2000) and Shim et al. (2009). These 
three Ha surveys are most similar to the likely observing mode 
of a JDEM/Euclid-like mission (slitless spectroscopy), and operate 
over a similar redshift range that will be pertinent to cosmology 
surveys. The adopted normalisation is 4>* — 1.37 x 10~' 5 Mpc -3 , 
and in Figure 1 we show how this compares to a range of observed 
LFs spanning the full redshift range < z < 2. Down to the lu- 
minosity corresponding to the flux limit likely to be practical in 
cosmology surveys (~10~ 16 ergs -1 cm~ 2 ), the simple model can 
replicate the observed space density of Ha emitters over 8 Gyr of 
cosmic time. At fainter limits, the uncertainty in the steepness of the 
faint-end slope will introduce further uncertainties that we ignore 
here, although we consider the effect of an evolving (steepening) 
faint-end in §2.2.1. 

Figure 2 shows another comparison to data that is more rele- 
vant for predictions for dark energy surveys - i.e. the redshift inte- 
grated counts as a function of limiting flux over 0.75 < z < 1.90 
(i.e. accessible in the near-IR). We compare the integrated counts 
derived from the Yan et al. (1999), Hopkins et al. (2000) and Shim 
et al. (2009) luminosity functions and the model. Note however, 
that these slitless surveys cover much smaller (<1 D°) areas than 
will be achievable with dedicated survey telescopes, and so suffer 
significantly from cosmic variance scatter - this could account for 
the scatter in the observations, and highlights the problem of the 
choice of normalisation mentioned above. The error-band in our 
model does not include the systematic uncertainty due to choice of 
4>*, but coincidentally spans the range of counts derived from the 
surveys shown in Figure 2. 

For convenience, we tabulate the predicted redshift distribu- 
tion dN/dz for a range of limiting fluxes, including uncertainties 
in Table 2. The distributions are plotted in Figure 3. Given the large 
scatter in the measured space density of Ha emitters determined 
from different surveys (see Hopkins et al. 2004 for a compilation), 
our adopted normalisation should be considered the best estimate 
'average'. However, when considering the feasibility of redshift 
surveys, the reader might want to adopt a more conservative es- 
timate of the density normalisation. If necessary, the reader can 
re-scale the predicted counts given in Table 2. We suggest that an 
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appropriate conservative lower limit to the counts could be taken as 
<f>* = 1 x 10" 3 Mpc- 3 .In§3 we discuss how the range of adopted 
normalisations affects our assessment of the feasibility of redshift 
surveys that aim to make cosmological measurements, and in the 
following section, we address further caveats that the reader should 
be aware of when applying this model. 

2.2 Caveats 

2.2.1 Evolution of the faint end slope 

Our model assumes a non-evolving faint end slope, with a — 
— 1.35 determined from local measurements (Gallego et al. 1995; 
Shioya et al. 2008). Both Yan et al. (1999) and Geach et al. 
(2008) fix this value of a in their fits of the luminosity function; 
the observations did not probe deep enough to constrain it. How- 
ever, there are hints that the relative abundance of galaxies with 
L < L* might increase towards early epochs (Reddy et al. 2008), 
with a as large as —1.6 at z ~ 2. What would be the ramifi- 
cations of a monotonically evolving (steepening) faint end slope 
out to z ~ 2 on our predicted counts? In Figure 3 we com- 
pare the redshift distributions for the fixed a model and the same 
model with a(z). We ignore the correlation between L*, <j>* and 
a for this analysis. At /n m > 10 -16 erg s _1 cm -2 , the counts 
(per redshift interval) predicted from the fixed a model are never 
less than ~85% of those derived from a steepening a model. At 
/iim > 5 x 10 -16 erg s _1 cm -2 the counts differ by only ~5%. 
This difference is smaller than the uncertainty on dN/dz, and so 
small enough to be ignored in this study. Needless to say, as high-z 
Ha studies probe deeper, past L* and can improve the constraint 
on a(z), the simple empirical model presented here could be re- 
vised accordingly. Finally, note that Hopkins et al. (2000) derive 
a faint-end slope of a — —1.6, which accounts for the turn-up in 
the integrated counts at / < 10~ 16 ergs -1 cm -2 (Fig. 2). The em- 
pirical model is not significantly different from the Hopkins et al. 
(2000) counts at brighter limits; re-enforcing that our assumption 
of a constant (local) faint-end slope is a reasonable baseline in this 
regime. 

2.2.2 Equivalent width cut 

An important feature of emission line surveys, and particularly nar- 
rowband surveys, is the inclusion of an equivalent width (EW) cut 
in the selection. Clearly this is an issue of sensitivity: galaxies with 
small EW are harder to detect and obtain reliable redshifts for. So 
naturally, dark energy surveys targeting emission-lines are biased 
towards galaxies with high equivalent widths, and against weak- 
emission lines and/or massive galaxies. In the model presented 
here we have assumed a fairly low EW cut, 10 A in the rest-frame. 
This cut will not significantly affect the predicted counts in the flux 
regime of interest. For example, according to the model of Baugh 
et al. (2005), at a flux limit fu m = 10~ 16 erg s _1 cm~ 2 , increasing 
the rest-frame equivalent width cut from 10A to 50A results in a 
drop in the number counts (integrated over 0.75 < z < 1.90 as in 
Fig. 2) of ~2%; the deficit is negligible at brighter limits. In prac- 
tice, redshift surveys will probably enforce an observed-frame cut 
of ~100A. 

Finally we note that the clustering properties of bright Ha 
emitters will be different from that of Ha emitters with low EW, 
or simply continuum- (e.g. JJ-band) selected galaxies. The latter 
should be more highly biased tracers of the mass distribution (see 
Orsi et al. 2009 in prep). 
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Figure 4. Prediction of contamination from [O II] emitters over the ob- 
served wavelength range 1-2/^m. To estimate the number of [O II] emitters, 
we have adapted the Ho count model, extrapolating to z = 4.4 (the redshift 
of [O II] at 2(j,m) making the assumption that all Ha emitters are also [O II] 
emitters, and these galaxies have a constant flux ratio of [O II]/Ha=0.62 
(Mouhcine et al. 2005). The contamination, expressed as a fraction of the 
total number of emitters detected, ranges from 1-13% in the range of limit- 
ing fluxes of practical interest. 

2.2.3 Contamination 

Emission-line surveys (aiming to detect a specific line; in this case 
Ha) are susceptible to contamination from galaxies with any strong 
emission lines at redshifts placing them in spectral range of the de- 
tector. At high-redshift this can be significantly problematic - for 
example, nearly two thirds of the potential z = 2.23 Ha emit- 
ters of Geach et al. (2008) selected with a narrowband at 2. 121/im 
were eliminated as low-redshift contaminants (e.g. Paa [z = 0.13], 
Pa/3 [z = 0.67], Fell [z = 0.3]). Higher redshift [O III]A5007 
can also contribute to the contamination. Geach et al. (2008) used 
further broad-band colour- and luminosity selections to select the 
2 = 2.23 candidates. Although most planned dark energy surveys 
will employ spectroscopy, one must still consider the potential for 
mis-identification of the Ha line in the large redshift ranges these 
surveys will probe. 

One could use the Ha model presented here to estimate the 
potential level of mis-identification of emission lines in spectral 
ranges likely to be employed in a slitless survey. For example, con- 
sider contamination from [O II] emitters at a rest-frame wavelength 
of 3727A. For a survey operating at 1-2/im, this means contam- 
ination from galaxies in the redshift range 1.7 < z < 4.4. If 
we assume that every Ha emitter is also an [O II] emitter, then 
we can estimate the expected number of objects in addition to the 
Ha emitters detected, assuming an attenuation due to the flux ra- 
tio [O Il]/Ha<l and intrinsic extinction -Afoul- In this example, 
we assume [O Il]/Ha=0.62 (measured from the 2dF Galaxy Red- 
shift survey, at z ~ 0.06; Mouhcine et al. 2005). Note that this 
ratio has not been corrected for the relative intrinsic extinction, and 
so this prediction should reflect the actual number of galaxies a 
flux-limited survey can expect to detecQ. As a fraction of the total 
number of emitters detected, the contamination from [O II] emitters 

1 Although strictly our model for the abundance of Ha emitters only ex- 
tends to 2 ~ 2, we assume the fixed evolution extends to 2 = 4.4. If the 
number of Ha emitters is actually gently declining at z > 2, this contami- 
nation estimate should be considered a conservative upper limit. 
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Table 2. Redshift distributions dN/dz (per square degree, calculated in bins of width Sz = 0.1, centred on the value given in the first column) for a range of 
limiting fluxes derived from the empirical model (also see Figure 3). For reference, we provide the range of Galactic extinctions at the observed wavelength 
of Ha, derived from the maps of Schlegel, Finkbeiner & Davis (1998). The predicted counts include intrinsic extinction in the Ha emitters, but the Galactic 
reddening will vary as a function of sky position. Although this has a negligible (few per cent) impact on the model dN/dz, we include it here as a guide. 
The counts listed here are calculated for a space density normalisation of <j>* = 1.37 X 10~ 3 Mpc~ 3 which is the 'average' space density of Ha emitters 
determined by several slitless surveys at z ~ 1 - similar to the Euclid and JDEM satellite survey concepts. The reader can re-scale these counts to alternative 
normalisations if desired: for a more conservative estimate of the counts, we recommend a lower density normalisation <f>* = 1 X 10 — 3 Mpc - 3 , however as 
we show in §3, this choice does not have a significant impact on the predicted power of a galaxy redshift survey. 
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ranges between 13% for a limiting flux of 1CF 16 ergs -1 cm -2 , to 
~1% for 5 x 10~ 16 ergs -1 cm -2 . A plot of the decline in con- 
tamination as a function of limiting flux is shown for reference in 
Figure 4. 

There are two simple ways to mitigate contamination. Perhaps 
the most efficient way to identify Ha is to resolve the [N II]A6583 
line (offset AA = 20A from Ha). Identifying this pair of lines 
is a useful discriminant between Ha and 'contaminant' lines, and 
so dark energy surveys should aim for a spectral resolution of 
R > 500 to achieve this. Another aid to redshift determination 
is the new generation of all sky ground based photometric surveys 
(e.g. PanSTARRS, Large Synoptic Survey Telescope). These sur- 
veys will provide optical photometry of many of the sources de- 
tected in the dark energy surveys; in conjunction with the near- 
IR photometry this will improve redshift estimates with a photo- z 
technique. 



2.2.4 Extinction 

The high redshift Ha surveys described in this work have not 
been corrected for intrinsic dust extinction, although when deriv- 
ing star formation rates, many authors tend to apply a canonical 
Aua = 1 mag unless some better estimate exists. The predicted 
number counts in our simple model include this intrinsic extinc- 
tion, such that if the extinction properties of the Ha emitters in 



the surveys described in Table 1 are relatively constant over a wide 
range of redshift, then the predicted counts can be taken as a re- 
liable representation of the expected yield even considering inter- 
nal extinction. However, all sky surveys (even ones that exclude 
the Galactic plane) will encounter a range of foreground Galactic 
extinction. Despite Ha being redshifted into the near-infrared at 
2 > 0.5, where reddening is fairly negligible, for completeness we 
consider here whether this could impact the predicted counts. 

Taking the all-sky dust maps of Schlegel, Finkbeiner & Davis 
(1998fl we evaluate the V-band extinction for Galactic latitudes 
\b\ > 20°, and extrapolate this to the observed wavelength of Ha, 
A = (1 + z) x 6563A out to z = 2.2 assuming a Rv = 3.1 red- 
dening law for the Galaxy (Cardelli et al. 1989; O'Donnell 1994). 
For reference, we summarise the average and range of reddenings 
for each redshift bin in Table 2. Of course, at longer wavelengths 
(in other words, Ha observed at higher-redshifts) reddening has an 
ever decreasing impact on the effective flux limit: at z > 0.5 the 
maximum Ah is never more than 0.2 mag, and the average is al- 
ways <0.03 mag. 

Since the regions of 'high' reddening represent a small frac- 
tion of the extragalactic sky, Galactic reddening has a minor 
(though redshift dependent) impact on the predicted counts. For 
example, modelling the variation in Ana over the full |6| > 20° 

2 irsa.ipac.caltech. edu/ applications /DUST/ 
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Figure 5. Predictions for the effective power of a galaxy redshift sur- 
vey, expressed in terms of the shot-noise parameter nP evaluated at 
k = 0.2/iMpc -1 (approximately the peak of the BAO signal). Fixed- 
time redshift surveys should aim for the sweet-spot of nPo.2 = 1 to ob- 
tain maximum power from the survey. We show the predicted np).2 for 
limiting fluxes of l-5x 10 — 16 erg s — 1 cm - 2 , and three survey 'efficien- 
cies' (e: the actual sampling of the Ha population due to the success rate 
of the survey). Note the clear degeneracy between survey efficiency and 
flux limit. The solid lines show the predictions for our 'average' model 
tf>* normalisation, but we also show the predicted nPo.2 for a more con- 
servative normalisation, </>* = 10 -3 Mpc -3 (for clarity only shown for 
/lim = 10 -16 ergs -1 cm -2 ). The conclusion to draw from this plot is 
that Ha surveys should be aiming for flux limits of ~10 — 16 erg s — 1 cm -2 ; 
beyond z ~ 1 the redshift yield goes into sharp decline, with severe conse- 
quences for nPo.2- 

sky, at z = 0.5 there is only a 2% decline in diV/dz; a smaller 
variation than the uncertainty of our model - we ignore its effects. 



3 IMPLICATIONS FOR REDSHIFT SURVEYS 

Dark energy surveys that aim to detect BAOs and measure red- 
shift distortions in galaxy clustering could target Ha emitters in 
all-sky near-infrared surveys, most likely utilising grisms for slit- 
less spectroscopy (e.g. McCarthy et al. 1999, Fig. 2). The key 
issue for these surveys is the ability to measure sufficient num- 
bers of redshifts for an accurate assessment of w(z) and f(z). 
Let us consider a hypothetical example: a slitless survey from a 
space platform with a wavelength coverage of l-2^im, and a spec- 
tral resolution of R > 500. This range gives access to Ha at 
0.5 < z < 2, with sufficient resolution to resolve [N II]A6583 
at /iim > 10~ 16 erg s -1 cm -2 . Aside from the slight modifica- 
tion to nominal limiting flux due to Galactic extinction (§2.2.4), 
there should be an additional modification to predicted counts due 
to some non-unity efficiency factor e (the ratio of the number of 
successfully measured redshifts, to the total number of measurable 



redshifts at a given flux limit). This will inevitably vary as a func- 
tion of flux, equivalent width, and so on). Including some assump- 
tion for e, how optimistic can we be about measurements of w(z) 
and f(z) in redshift surveys? 

A precise measurement of w(z) or f(z) requires an accurate 
measurement of the power spectrum, P{k). The uncertainty with 
which P(k) can be measured from a given galaxy survey depends 
on the number density of galaxies and the volume of the survey. If 
the number density is low, then the errors are dominated by shot 
noise. If it is high, then cosmic variance (i.e. the volume of the sur- 
vey) dominates the error budget. To see this, note that the effective 
volume of a survey is given by Feldman, Kaiser & Peacock (1994) 
as: 



/dV 


n(r)P 


1 + n(r)P 





(3) 



where n(r) is the comoving number density of the sample at loca- 
tion r. For small n, V a s oc n and the signal is shot noise dominated. 
For large n, V g = V, where V is the physical volume of the sur- 
vey, which limits the signal. For a sample with a fixed total num- 
ber of galaxies iVg a i = nV, and for a power spectrum P, setting 
d\4ff/dV = requires nP = 1. In this situation we see that the 
effective volume reaches a maximum when nP = 1. This 'sweet- 
spot' is often used as a design aim for fixed integration-time and/or 
volume limited galaxy redshift surveys, with P0.2 = (P(fc)), cal- 
culated for k = 0.2 /i Mpc -1 . This scale is approximately the limit 
of the quasi-linear regime, and this also gives an indication of the 
strength of the clustering signal on the linear scales carrying the 
redshift-distortion information. 

Future surveys will often be limited by the extragalactic sky 
area they can observe. For surveys using a single ground-based 
telescope (such as BOSS) this is of order ~10 4 square degrees, 
while for a space based platform (such as Euclid or JDEM) or a 
survey using a pair of telescope in different hemispheres, this is 
~(2 — 3) x 10 4 square degrees. In this situation, the volume that 
can be surveyed in the interesting redshift range is limited, and the 
only way of gaining signal is to push to higher galaxy number den- 
sities. It is therefore important to consider values of nPo.2 > 1. 

In Figure 5 we show the predicted nPo.2 as a function 
of redshift, for a range of (nominal) limiting fluxes (1-5) x 
10~ 16 ergs -1 cm -2 . As well as the ideal case, with an efficiency 
factor e — 1 (that is, one correctly identifies all the Ha emitters 
above the survey flux limit in every pointing), we show the effect on 
nPo.2 for a 50% and 25% efficiency. Note that we have assumed a 
model for the luminosity-dependent evolution of bias for Ha emit- 
ters from Orsi et al. (2009) such that P ga i = Pr>Mb(z, Lua) 2 - 
The Ha population is generated by the semi-analytic prescription 
GALFORM (Baugh et al. 2005). Since in the semi-analytic model 
one can ask what dark matter halo hosts a given galaxy, Orsi et 
al. estimate the galaxy bias for a given Ha luminosity by averag- 
ing over the halos that host selected Ha emitters. The model bias 
for Ha emitters at z ~ 2 agrees well with the value derived by 
Geach et al. (2008) from the projected two-point correlation func- 
tion. Note that we have applied the same rest-frame EW cut as ap- 
plied throughout this work, and interpolated the b(z, Lhc) as nec- 
essary. As a guide, the range of bias applied over < z < 2 for the 
luminosities corresponding to the limiting fluxes considered here is 
0-9 < b < 1.7. 

Obviously one would always strive for maximum efficiency 
and depth, but this is not a practical possibility: there will always 
be redshift attrition resulting in e < 1. This inefficiency has the 
same impact as increasing the effective limiting flux of the survey. 
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Losing counts has a serious impact on the survey power; even at the 
faintest limit likely to be practicable, 10 -16 erg s _1 cm -2 , a 'per- 
fect' survey struggles to achieve nPo.2 = 1 at z = 2. Assuming 
the more likely case of e = 0.5, one can comfortably achieve the 
required nPo.2 out to z — 1, even with fairly conservative flux lim- 
its. At higher redshifts this becomes increasingly observationally 
expensive. Re-visiting the caveat of model normalisation described 
in §2.1, on Figure 5 we also show the more conservative case the 
reader might choose to adopt. Obviously a shift in normalisation 
simply translates the predicted nPo.2 up or down. It is worth noting 
that the conservative counts are within the la band of uncertainty 
of the average model normalisation at z ^ 1, and so our conclu- 
sions about the power of redshift surveys as a function of limiting 
flux and efficiency are unchanged. 

One way to boost performance would be to employ Digital 
Micro-mirror Devices (DMDs), rather than traditional slitless spec- 
troscopy. For a fixed telescope diameter and integration time, with 
DMD-slit spectroscopy one reaches ~2.5 mag deeper in the contin- 
uum, due to the strong reduction of the sky background compared 
to slitless spectroscopy. This allows the detection of several spectral 
features in each spectrum (absorption and emission lines) and the 
consequent identification of all galaxy types (early-type and star- 
forming systems). Moreover, thanks to improved sensitivity and 
the lack of the 'spectral confusion' problem due to the overlap of 
spectra of different objects (the traditional Achille's heel of slit- 
less spectroscopy), the redshift success rate e is much higher (up to 
>90%, see Cimatti et al. 2009). 



4 COSMOLOGICAL NEAR-IR SURVEYS FROM THE 
GROUND 

To be competitive with space-platforms targeting Ho emitters at 
z > 0.5, ground-based near-IR BAO surveys should also be aim- 
ing for limiting fluxes of 10~ 16 erg s _1 cm -2 , but there are extra 
observational challenges - not least the deleterious effect of the 
atmosphere in the near-IR. Approximately 30% of the 1-2/im win- 
dow has an atmospheric transmission of <80%, mainly affecting 
Ha in the redshift ranges 1 < z < 1.3 and 1.7 < z < 2.1. In 
addition, near-IR observations from the ground must also contend 
with forest of OH-airglow: even at R ~ 2000 less than half of the 
near-IR spectral range is free from OH line emission, although new 
OH suppression technologies could partly mitigate this effect. 

On the basis of areal coverage, ground-based near-IR BAO 
survey will never be competitive with a Euclid/JDEM-like mission. 
Modern wide-field near-IR spectrographs deploy fibres on individ- 
ual targets, and this presents a significant disadvantage compared to 
the slitless approach of the space missions: one must select targets 
prior to observation (in some sense the problem is reversed in the 
slitless case). Typically this will require the target fields to be com- 
plemented by multi-colour broad-band photometry, deep enough to 
provide an estimate of redshift. Note that the consequence for mis- 
identifying targets is a strong hit to the efficiency parameter e. 

In the event of a dedicated space-based near-infrared dark en- 
ergy survey going ahead, one could argue that a more efficient use 
of ground based multi-object spectrographs in the near-IR would be 
to complement the wider cosmological surveys by providing more 
detailed follow-up observations of a sub-sample of line-emitters. 
This has the advantage of side-stepping the issue of target selec- 
tion, since the sample would already be 'sanitized' by the cosmol- 
ogy survey. Such a symbiosis between space and ground would be 
an efficient use of resources since: (a) the ground facilities would 



target known line-emitters, and therefore rapidly build up a large 
sample of spectroscopic observations for high-redshift galaxies in 
more detail than can be achieved from the space platforms; and (b) 
the complementary observations could help to better characterise 
contamination from other line emitters (as discussed in §2.2.3), thus 
feeding back information to the cosmological survey. Combining 
surveys in this way could serve to satisfy two groups of researchers: 
those interested in the astrophysics of galaxies at high redshift, and 
those concerned with cosmological measurements. 



5 SUMMARY AND FINAL REMARKS 

We have presented a simple prescription for the prediction of the 
abundance of Ha emitters over < z < 2, based on empirical 
data. The model is simplistic, due to limited available data; it as- 
sumes a fixed space density, fixed faint end slope, and only L* evo- 
lution out to z = 1.3. There is no luminosity evolution to higher 
redshifts, consistent with current Ha observations at this redshift. 
Despite its simplicity, the model adequately mimics the observed 
luminosity functions of a range of Ha surveys (including a mix- 
ture of spectroscopic, grism and narrowband strategies). Using the 
luminosity function model as a basis, we predict the redshift dis- 
tribution of Ha emitters corresponding to a spectral coverage that 
extends to 2/iin. 

Our results have particular relevance to dark energy experi- 
ments attempting to measure cosmological information from the 
power spectrum of galaxies detected in all-sky Ha surveys in the 
near-IR. We use the parameter nPo.2 as a measure of the effective- 
ness of a redshift survey, and make predictions for this value for 
a range of redshift, limiting flux and success rate (i.e. efficiency). 
To achieve nPo.2 = 1 out to z — 2, emission-line surveys should 
be aiming for limiting fluxes of ~10 -16 erg s _1 cm~ 2 . However, 
this estimate is reliant on a high success rate of the sampling of the 
Ha population: redshift surveys need to aim for high-efficiencies, 
since any decline in redshift yield (i.e. failing to obtain redshifts 
for detections) has the same effect on nPo.2 as increasing the flux 
limit (illustrated in Figure 5 of this work). Assuming a more likely 
situation of 50% efficiency, a realistic target for proposed surveys 
is nPo.2 = 1 at z = 1.5. At higher redshifts the sharply declin- 
ing number counts have a severe effect on one's ability to measure 
w(z) at the desired precision. 
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